Quantum WDM fermions and gravitation determine the observed galaxy structures 
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Quantum mechanics is necessary to compute galaxy structures at kpc scales and below. This is so 
because near the galaxy center, at scales below 10 — 100 pc, warm dark matter (WDM) quantum 
effects are important: observations show that the interparticle distance is of the order of, or smaller 
than the de Broglie wavelength for WDM. This explains why all classical (non-quantum) WDM 
iV-body simulations fail to explain galactic cores and their sizes. We describe fermionic WDM 
galaxies in an analytic semiclassical framework based on the Thomas-Fermi approach, we resolve it 
numerically and find the main physical galaxy magnitudes: mass, halo radius, phase-space density, 
velocity dispersion, fully consistent with observations, including compact dwarf galaxies. Namely, 
fermionic WDM treated quantum mechanically, as it must be, reproduces the observed galaxy DM 
cores and their sizes. [In addition, as is known, WDM simulations produce the right DM structures 
in agreement with observations for scales > kpc] . We show that compact dwarf galaxies are natural 
quantum macroscopic objects supported against gravity by the fermionic WDM quantum pressure 
(quantum degenerate fermions) with a minimal galaxy mass and minimal velocity dispersion. 
Interestingly enough, the minimal galaxy mass implies a minimal mass m m in for the WDM particle. 
The lightest known dwarf galaxy (Willman I) implies m > m m in = 1.91 keV. These results and the 
observed halo radius and mass of the compact galaxies provide further indication that the WDM 
particle mass m is approximately around 2 keV. 

PACS numbers: 95.35.+d, 98.52.-b, 98.56.Wm, 98.62. Gq 
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I. INTRODUCTION AND SUMMARY OF RESULTS 



Dark matter (DM) is the main component of galaxies, especially of dwarf galaxies which are almost exclusively 
formed by DM. It thus appears that the study of galaxy properties is an excellent way to disentangle the nature of 
DM. 
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Cold DM (CDM) produces an overabundance of substructures below the ~ 50 kpc till very small scales ~ 0.005 
pc which constitutes, as is well known, one of the most serious drawbacks for CDM. On the contrary, warm DM 
(WDM), that is, DM particles with mass in the keV scale, produces DM structures in the range of scales < 50 kpc 
in agreement with observations. In WDM structure formation, substructures below the free-streaming scale ~ 50 
kpc are not formed, contrary to the case of CDM. This conclusion for WDM based on the linear theory is robustcly 
confirmed by A-body simulations by different groups For scales larger than 50 kpc, WDM yields the same results 
than CDM and agrees with all the observations: small scale as well as large scale structure observations and CMB 
anisotropy observations. 

Astronomical observations show that the DM galaxy density profiles are cored till scales below the kpc 2 4]. On 
the other hand, TV-body CDM simulations exhibit cusped density profiles with a typical 1/r behaviour near the ga laxy 
center r = 0. Classical TV-body WDM simulations exhibit cusps or small cores smaller than the observed cores [a.[la|. 

Numerical calculations based on the spherically symmetric Vlasov-Poisson equation from the Larson moment 
expansion , as well as on the exact dynamics of the associated N-body system, have confirmed these findings [|[ . 

A direct way to see whether a system of particles has a classical or quantum nature is to compare the particle de 
Broglie wavelength XdB with the inter-particle distance d. We investigate this issue in sec. II Al and express the ratio 
of the two lengths as 

n = — = h — , Q h = — . 



where Qh is the DM phase space density, ph and a being the halo density and the velocity dispersion, respectively. 
Values of 1Z much smaller than unity correspond to classical physics while 1Z ~ 1 or larger corresponds to the quantum 
regime. The observed values of Qh from Table U yields 1Z in the range 

2xl0- 3 <ft (-^-f <1.4. (1.1) 
VkeV / 

The larger value of 1Z is for ultracompact dwarf galaxies and the smaller value of 1Z is for large spirals. The values 
of 1Z around unity clearly imply (and solely from observations) that compact dwarf galaxies are natural macroscopic 
quantum objects for WDM. 

WDM fermions always provide a non-zero pressure of quantum nature. By balancing this quantum pressure with 
the gravitation pressure, we find theoretical values for the total mass M ~ 10 6 Mq, the radius R ~ 30 pc and the 
velocity dispersion a ~ 2 km/s which are consistent with the observations of compact dwarf galaxies (see Table HJ. 
These results back the idea that dwarf galaxies are supported by the fermionic WDM quantum pressure. 

Classical A-body simulations for DM are not valid at scales where the interparticle distance becomes of the order 
of, or smaller than, the de Broglie wavelength. This precisely happens near the galaxy center for WDM particles with 
mass m ~ keV. Therefore, the results of classical (non-quantum) WDM TV-body simulations at such scales are not 
valid. This explains why all classical WDM A-body simulations fail to correctly reproduce the observed galaxy cores 
Classical A-body simulations are reliable for CDM because the mass of the CDM particle is large enough (m ~ 
GeV) and then TZ -C 1. Such reliable CDM simulations always produce cusped galaxy profiles in contradiction with 
the observations. 

Richard P. Feynman foresaw the necessity to include quantum physics in simulations § : 

"I'm not happy with all the analyses that go with just the classical theory, because nature isn't classical, dammit, and 
if you want to make a simulation of nature, you 'd better make it quantum mechanical, and by golly it 's a wonderful 
problem, because it doesn't look so easy. " 

We treat here the self-gravitating fermionic DM in the Thomas-Fermi approximation. In this approach, the DM 
chemical potential /j,(r) = po — m 4>{r), where /j,q is a constant and 4>{r) the gravitational potential, obeys the 
selfconsistent Poisson equation 



d n 2 dfi _ 
dr 2 r dr it 2 h 3 



p 2 dp f (e(p) - n(r)) (1.2) 



(i 



where G is Newton's gravitational constant, p is the DM particle momentum, e(p) = p 2 /(2m) is the DM particle 
kinetic energy and f(E) is the energy distribution function. This is a semiclassical gravitational approach to determine 
selfconsistently the gravitational potential of the fermionic WDM given its distribution function /. 
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The boundary condition /i'(0) = at the origin guarantees bounded central DM mass densities. 

The distribution function f(E) (or a more general form with a dependence on other constants of the motion besides 
E) is determined by the DM evolution since decoupling. Such quantum dynamical calculation is beyond the scope 
of the present paper. f(E) can be modelized for instance by the equilibrium Fermi-Dirac distributions or by out of 
equilibrium distributions. The galaxy magnitudes turn to be rather insensible to whether one chooses equilibrium or 
out of equilibrium distributions f(E) Q. 

We get a one parameter family of solutions of eas. (|1.2l) parametrized by the value of the chemical potential at the 
origin /«(0) that can be expressed in terms of the phase-space density at the origin Q(0). Large positive values of /i(0) 
correspond to most compact object (fermions in the quantum degenerate limit), while large negative values of fj,(0) 
yield dilute objects (classical limit). 

We show that the Thomas-Fermi equation implies the local equation of state 

P(r) = cr 2 (r) p(r) and the hydrostatic equilibrium equation — — h p(r) -f- = . 

dr dr 

This local equation of state generalizes the local perfect fluid equation for r-dependent velocity v(r). 

The numerical resolution of eqs. (|1.2p for the whole range of the chemical potential at the origin /z(0) yields the 
physical galaxy magnitudes, such as mass, halo radius, phase-space density and velocity dispersion all fully compatible 
with observations including compact dwarf galaxies as can be seen from figs. [T] and [5] and Table HI 

Approaching the classical diluted limit yields larger and larger halo radii, galaxy masses and velocity dispersions. 
Their maximum values arc limited by the initial conditions provided by the primordial power spectrum which deter- 
mines the sizes and masses of the galaxies formed. The phase space density decreases from its maximum value for the 
compact dwarf galaxies corresponding to the degenerate fermions limit till its smallest value for large galaxies (spirals 
and ellipticals) corresponding to the classical dilute regime. The theoretical values for the core radius r^, the core 
galaxy mass Mh and the velocity dispersion er(0) obtained in the Thomas-Fermi framework vary very little with the 
specific form of the distribution function f(E). 

We display in fig. [T] the density and velocity profiles, p(r)/p(0) and a(r)/a(0) obtained in the Thomas- Fermi 
approach for different values of the chemical potential at the origin /i(0). Large positive values of /i(0) correspond to 
the compact galaxies, while negative values of /i(0) correspond to the classical regime describing spiral and elliptical 
galaxies. All density profiles are cored. The sizes of the cores 77, are in agreement with the observations, from the 
compact galaxies where r% ~ 35 pc till the spiral and elliptical galaxies where r/, ~ 0.2 — 60 kpc. The larger and 
positive is p.(0), the smaller is the core. The minimal one arises in the degenerate case p(0) — > +oo (compact dwarf 
galaxies) . 

In the left panel of fig. [5J we plot the (dimensionless) galaxy phase-space density ft 3 Q(0)/(keV) 4 obtained from 
the numerical resolution of the Thomas- Fermi ecis. (|1.2p for WDM fermions of mass m = 1 and 2 keV. The observed 
values of ft 3 Q^/(keV) 4 and rt from Table U are also depicted. The theoretical Thomas-Fermi curves in fig. [5] appear 
slightly below the observational data in all the range of galaxies most likely because the observed values Qh (derived 
from the stars' velocity dispersion) are in fact upper bounds for the DM Qh- 

In the right panel of fig. [2J we plot the galaxy masses (Mh/ Mq)^J M©/pc 3 po obtained from the numerical resolution 
of the Thomas- Fermi eqs. (|1.2p . where Mh is the mass inside the halo radius and po the observed central mass density. 
We observe a good agreement between the Thomas-Fermi results and the observations in all the range of galaxies for 
a DM particle mass m around 2 keV. Notice that the error bars of the observational data are not reported here but 
they are at least about 10 — 20%. 

For degenerate fermions [p(0) — > +oo] the halo radius, the velocity dispersion and the galaxy mass take their 
minimum values. These minimum values eqs. (|2.25p and (|2.25l) obtained in the Thomas- Fermi approach are similar 
to the estimates provided in sec. II Bl through a simple balancing argument between the gravitational and quantum 
pressures. 

The masses of compact dwarf galaxies dominated by DM must be larger than this minimum mass Mh, m in- The 
lightest known galaxy of this kind is Willman I (see Table|T|. Imposing Mh,min < Mwaiman i — 2.9 10 4 Mq provides 
a minimal mass (a lower bound) for the WDM particle: 

m > m m i n = 1.91 keV . (1.3) 

This minimal WDM mass is independent of the WDM particle physics model. m m i n is an universal value irrespective 
of the shape of the distribution function f(E) in the non-degenerate regime. 
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TABLE I: Observed values r^, a, y/Qh, p(0) and Mh covering from ultracompact galaxies to large spiral galaxies from 
refs. @, i, [H 51 [HHI| . The phase space density is larger for smaller galaxies, both in mass and size. Notice that the phase 
space density is obtained from the stars velocity dispersion which is ex pec ted to be smaller than the DM velocity dispersion. 
Therefore, the reported Qh are in fact upper bounds to the true values [lq ]. 

Interestingly enough, comparison of the theoretically derived galaxy masses with the galaxy data plotted in fig. 
[2] indicates a WDM particle mass m approximately around 2 keV in agreement with earlier estimations floL ll2T | . If 
the WDM particle mass would bera> 1 keV, an overabundance of small galaxies (small scale structures) without 
observable counterpart would appear. 

In summary, the theoretical Thomas-Fermi results are fully consistent with all the observations including dwarf 
compact galaxies as can be seen from figs. [T] and [2] It is highly remarkable that in the context of fermionic WDM 
the simple stationary quantum description provided by this semiclassical framework is able to reproduce such broad 
variety of galaxies. 

These results indicate that fermionic WDM treated quantum mechanically (even approximately) is fully consistent 
with the observed galaxy properties including the DM core sizes. That is, quantum physics and the associated 
quantum pressure, rule out galaxy cusps for fermionic WDM and provide the right sized observed cores. 

We have not considered baryons in the present analysis of the galaxies. This is fully justified for dwarf compact 
galaxies which are composed today 99.99% of dark matter [HI HH [Hj] . ^ n l ar S e galaxies the baryon fraction can reach 
values up to 1 - 3 % [l4[. We have also ignored supermassive central black holes which appear in large spiral galaxies. 
In any case, it must be noticed that the central black hole mass is at most ~ 1CU 3 of the mass of the bulge. 
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Fermionic WDM by itself produce galaxies and structures in agreement with observations. Therefore, the effect of 
including baryons is expected to be a correction to the pure WDM results, consistent with the fact that dark matter 
is in average six times more abundant than baryons. 

We use units such that the speed of light is c = 1 throughout this paper. 

A. Dwarf galaxies as WDM quantum macroscopic objects 

In order to determine whether a system of particles has a classical or quantum nature we should compare the 
particle de Broglie wavelength with the interparticle distance. 

The de Broglie wavelength of DM particles in a galaxy can be expressed as 

^dB = —— , (1.4) 

to a 

where a is the velocity dispersion, while the average interparticle distance d can be estimated as 

"=y ■ (i - 5) 

where ph is the average density in the galaxy core. We can measure the classical or quantum character of the system 
by considering the ratio 

d 

By using the phase-space density 

n - ph 

Qh = 

and eqs. (jl.4|) - (|1.5|> . TZ can be expressed as 

*-»($)'. (1-6) 

Notice that TZ as well as Qh are invariant under the expansion of the universe because both X^b and d scale with the 
expansion scale factor. TZ and Qh evolve by nonlinear gravitational relaxation. 
Using now the observed values of Qh from Table [J yields TZ in the range 

2 x 1(T 3 < U {^rj < 1.4 (1.7) 

The larger value of TZ is for ultracompact dwarfs while the smaller value of TZ is for big spirals. 

The ratio TZ around unity clearly implies a macroscopic quantum object. Notice that TZ expresses solely in terms of 
Q and hence (h 3 Q/m A ) measures how quantum or classical is the system, here, the galaxy. Therefore, we conclude 
solely from observations that compact dwarf galaxies are natural macroscopic quantum objects for WDM. 

B. Dwarf Galaxies supported by WDM fermionic quantum pressure 

For an order-of-magnitude estimate, let us consider a halo of mass M and radius R of fermionic matter. Each 
fermion can be considered inside a cell of size Ax ~ l/ns and therefore has a momentum 

h * i 

P n*> — ~ h 77,3 . 

Ax 

The associated quantum pressure P q (flux of the momentum) has the value 

P„=nap~h<rn3 = — n 3 . (1.8) 

TO 
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where a is the mean velocity given by 



We estimate the number density as 



and we use that p = m a to obtain from eq. 



p h i 
- = -n' 

m m 



M 
1 7r R 3 m 



ft? f 3M V 

q ~^R? \Jr7^) • (1 ' 9) 

On the other hand, as is well known, galaxy formation as all structure formation in the Universe is driven by 
gravitational physics. The system will be in dynamical equilibrium if this quantum pressure is balanced by the 
gravitational pressure 

G M 2 1 

Pa = gravitational force/area = - — x (1.10) 

5 ' R 2 At: R 2 V ' 

Equating P q = Pq from eqs. (|1.10[) - (|1.9[) yields the following expressions for the size R and the velocity a in terms of 
the mass M of the system and the mass m of the particles,: 



3S h 2 /10 b M, 



/kevy 



R = — — rr — — 8—ttt = 10 -61 ■ ■ ■ pc — T-r-^ ) , (1.11 

(4 7r)3 G ma Ms V M J \ m J 

1 4 / 2 

ATT Y G 4 , ,2 _ _„ km / m \ s ( M 



, _ ma M3 = 11.62... ' — . (1.12) 

81 J h s VkcV/ V 10 6 M eJ 

Notice that the values of M, R and a are consistent with the observed values of dwarf galaxies. Namely, for M of 
the order 10 6 M Q (which is a typical mass value for dwarf galaxies), R and a give the correct order of magnitude for 
the size and velocity dispersion of dwarf galaxies (see Table |T| for a WDM particles mass in the keV scale. 

These results back the idea that dwarf galaxies are supported by the fermionic WDM quantum pressure eq. (|1.9[) . 

It is useful to express the above quantities in terms of the density p, as follows 



M _ •* JjZ _ 0.7075 . . M S I (i£X)\ « . « ^ _ 25 . 66 . . . pc (_««.)■ (i£X Y , 

2 777 4 V I" G J y Mq \ m J 2 V 7T G 777 3 pe \PP C / V m J 

t / p \s km / pc 3 V /keVV t2 M Q / pc 3 V /keV\^ 

& = h ( ——r ) = 1.988 ... — p^- ,p=ft 2 il^ = 0.04399 . . . — % [p ¥—) . 1.13 

VmV s V M Q ) \m J q 777 f kpc 3 \ Mq J \mj v ' 

R and M are typical semiclassical gravitational quantities involving both G and h. The particle velocity a and the 
pressure P q are of purely quantum mechanical origin. 



II. GENERAL GALAXY PROPERTIES FROM QUANTUM FERMIONIC WDM IN THE 

THOMAS-FERMI APPROACH 

We consider a single DM halo in the late stages of structure formation when DM particles composing it are non- 
rclativistic and their phase-space distribution function f(t, r,p) is relaxing to a time-independent form, at least for 
r not too far from the halo center. In the Thomas-Fermi approach such a time-independent form is taken to be a 
energy distribution function f(E) of the conserved single-particle energy E = p 2 /(2m) — p, where m is the mass of 
the DM particle and p, is the chemical potential 

/i(r) = po -m(f)(r) (2.1) 
with <^(r) the gravitational potential and p some constant. 
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We consider the spherical symmetric case where the Poisson equation for (j)(r) takes the form 

d 2 a 2 da , . . . 

^ + -£ = -4«Gmp(r), (2.2) 

where G is Newton's constant and p(r) is the DM mass density. In turn, p(r) is expressed here as a function of p{r) 
through the standard integral of the DM phase-space distribution function over the momentum 

where g is the number of internal degrees of freedom of the DM particle, with g = 1 for Majorana fermions and g = 2 
for Dirac fermions. For defmiteness, we will take g — 2 in the sequel. 

Another standard integral of the DM phase-space distribution function is the pressure 
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and from p(r) and P(r) other quantities of interest, such as the velocity dispersion cr(r) and the phase-space density 
Q{r) can be determined as 

-»=^4 , qm=4V (2.5) 

p[r) a a \r) 

We see that p(r) fully characterizes the fermionic DM halo in this Thomas-Fermi framework. The chemical potential 
is monotonically decreasing in r since eq. (|2.2[) implies 

^ = _GmM(r) M (r) = An T dr' r' 2 p(r') . (2.6) 

dr r z J Q 

Moreover, the fermionic DM mass density p is bounded at the origin due to the Pauli principle Q, and therefore the 
proper boundary condition at the origin is 

f(0)=0. (2.7) 

Also notice from eq. (|2.ip that, by assuming 0(0) = 0, then po = u(0), which plays the role of a Lagrange multiplier 
to be eventually eliminated in favor of po = p(0), the mass density at the origin. 

Eqs. (|2.2[ ) and (|2.3I) provide an ordinary nonlinear differential equation that determines selfconsistently the chemical 
potential p(r) and constitutes the Thomas-Fermi semi-classical approach. Fermionic DM in this approach has been 
previously considered in ref. [l3j . 

In this semi-classical framework the stationary energy distribution function f(E) must be assigned beforehand. In 
a full-fledged treatment one would rather solve the cosmological DM evolution since decoupling till today, including 
the quantum dynamical effects which become important in the non-linear stage and close enough to the origin. Such 
a quantum dynamical calculation is beyond the scope of the present paper. 

From eq. (|2.4p and (|2.5p we derive the local equation of state: 

P(r)=a 2 (r) p(r) , (2.8) 

and the hydrostatic equilibrium equation 

dP , . dd> ,„ „, 

-7- + pM-f=0. (2.9) 
dr dr 

The local equation of state eq. (|2.8p generalizes the local perfect fluid equation of state for r-dependent velocity v(r). 
As we see below, the perfect fluid equation of state is recovered both in the classical dilute limit and in the quantum 
degenerate limit. 

Eliminating P(r) between eqs. (|2.8p and (|2 .9[) and integrating on r gives 

f r dr' d(j) 
P( r ) = g2 ( ) ~ J ^V) ^7 
p(0) a 2 (r) 
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For constant v(r) this relation reduces to the baryotropic equation. Inserting this expression for p(r) in the Poisson's 
equation yields 

dr' d<j) 



d 2 6 2 d6 . „ cr 2 (0) / „2( r i\ a t i 
-r-7T + - -p = 4tt G m po -^-f e Jo <r (r j 
dr z r dr a z (r) 

This nonlinear equation generalizes the corresponding equation in the self-gravitating Boltzmann gas when cr 2 (r) is 
constant. 



We integrate the Thomas-Fermi nonlinear differential equations (|2.2[) - (|2.3|) from r = till the boundary r = R = 
i?200 ~ Rvir defined as the radius where the mass density equals 200 times the mean DM density. It is useful to define 
dimensionless variables £, £r, and the dimensionless one-parameter distribution function \E' as 

r = /o£ , R = hin , ii(r) = Eq u(£) , f(E)=*(J^j, (2.10) 

where Po is the characteristic one-particle energy of the DM halo and Iq is the characteristic length that emerges 
from the dynamical equations (|2.2j) - (|2.3|) : 



9tt V _ /keVV / pc 3 v " 



Z = —= — — -i? Pofr ' = 18.709 ... pc , (2.11) 

V8G \ m PoJ \ m J \ M oJ 

The self-consistent Thomas-Fermi equation (|2.2l) - (|2.3j) for takes the form 
d 2 v 2 dv _ I 2 {v) 



where 



, pfa) = 200 pdm , i/(0) = , vb = i/(0) , (2.12) 



J n (i/) = (n+l) / y n dy*(y'-v) , n=l,2,... , (2.13) 



and we use the integration variable y = p/y/2 m Eq. 

We find the main physical galaxy magnitudes, such as the mass density p(r), the average velocity of the particles 
(r) and the pressure P(r) (which are all r-dependent) as: 



v 



M r ) = Po r ; s , <r(r) = a 0) ■— — , P r = P 0) , . , 2.14 



m 4 /2P ^ 3/2 



(-Y 


Po 


V TO J 





5/3 

^4(^0) 



rn s „ frfo)] 1 ' 2 AeV\ 4 / 3 / pc 3 V' a km 



, 4/3 / pc^ 1; 

[/2K)] 5 / 6 I m J M y s 

As a consequence, the total mass Mr enclosed in the sphere of radius R and the phase space density Q(r) turn to be 
Mr = 4tt /^r 2 p(r) = An f/\ f* dxx 2 I 2 (u(x)) = 4tt f f & = 

JO M^o) Jo [M*W /3 



£r k'(6t)l /keV\ 4 / pc 3 4tt /P \ 3 , 

= Mo fekP (— ) r ' Mo " M0 U) = °- 82296 ■ • ■ x 10 ' (2 - 16) 

Q(r) = Q(0) 4^ 4^ , Q(0) = ^ m 4 ^ , VQ = m , po = p(0 ) . (2.17) 

We have systematically eliminated the energy scale Eq in terms of the central density pq. [We may also choose the 
density at another point r ^ 0]. Notice that Q(r) turns to be independent of Eq and therefore from pq. 
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Besides the virial galaxy radius R = Iq £r, one can define the core size of the halo by analogy with the Burkert 
density profile as 



P{ru) = 1 
Po 4 



r h =lo£h. (2-18) 



To explicitly solve eq. (12.12)) we need to specify the distribution function *S?(E/Eq). But many important properties 
of the Thomas-Fermi semi-classical approximation do not depend on the detailed form of the distribution function 
H?(E/E ). Indeed, a generic feature of a physically sensible one-parameter form &(E/E ) is that it should describe 
degenerate fermions for E —> 0. That is, fy(E/E ) should behave as the step function 8(—E) in such limit. In the 
opposite limit, ^(E/Eq) describes classical particles for (i/Eq — oo. As an example of distribution function, we 
consider the Fermi-Dirac distribution 

*ME/Eo) = ■ (2-19) 

The choice of ^fd may be justified near the origin, where relaxation to thermal equilibrium is conceivable. Far from 
the origin however, the Fermi-Dirac distribution as its classical counterpart, the isothermal sphere, produces a mass 
density tail 1/r 2 that overestimates the observed tails of the galaxy mass densities. Notice indeed that the classical 
regime (j,/Eq — > — oo may be attained for large distances r since eq. (|2.6l) indicates that /x(r) is always monotonically 
decreasing with r. 

More precisely, large positive values of the chemical potential at the origin vq ^> 1 correspond to the degenerate 
fermions limit which is the extreme quantum case and oppositely, vq <C — 1 gives the diluted limit which is the classical 
limit. In this classical limit the Thomas-Fermi equations (12.2l) - (|2.3p become the equations for a self-gravitating 
Boltzmann gas. 

We display in fig. [T]the density and velocity profiles. Namely, we plot p(r)/p and a(r)/a(0) as functions of r/R 
for vq = i>(0) — —5, 0, 5, 15, 25 and in the degenerate fermion limit vq — > +oo. The obtained fermion profiles are 
always cored. The sizes of the cores defined by eq. (|2.18[) are in agreement with the observations, from the compact 
galaxies where ~ 35 pc till the spiral and elliptical galaxies where ~ 0.2 — 60 kpc. The larger and positive is Uq, 
the smaller is the core. The minimal core size arises in the degenerate case vq — > +oo (compact dwarf galaxies). 

Eqs. (|2. 14|) - (|2.16|) cover the full range of physical galaxy situations from the quantum degenerate fermions (dwarf 
galaxies) until the dilute classical limit (spiral and elliptic galaxies) as we discuss in subsection III CI In addition, 
the galaxy length scale Zo and the galaxy mass Mr emerging from the Thomas-Fermi approach turn to be of the 
same order of magnitude of the Jeans' length and Jeans' mass, respectively which are here of semi-classical nature, 
containing h and G @. 



A. Quantum degenerate limit: compact dwarf galaxies and minimal galaxy mass 

The quantum degenerate limit is very instructive and we consider it now. In this limit » 1, the distribution 
function ^(E/E ) becomes a step function 6(—E) and all fermion states get filled till the Fermi level. Notice that 
the degenerate limit of the distribution function is always a step function 6{—E) irrespective of the shape of the 
non-degenerate fermion distribution function ^(E/Eq). Then, from eq. (l2.13[) 

In{v)deg = ^ , (2.20) 

The Thomas-Fermi equation (|2 . 1 2[) becomes 

^ + |§ = -^ 3/2 , m = »?(0) = 1, 1/(0) =0, (2.21) 

which we solve numerically for < £ < with £^ = 3.6537446 . . . the first zero of ry(£) and the dimensionless core 
radius £/, = 2.269587 . . .. One then finds for the mass density and total mass in the degenerate limit 

o3/2 T „5/2 /-ffl 
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M degt h = 1.856 . . . M (^pj y po |£ , (2.22) 

where the numerical result £fj |?/(£i?)l = 2.71405512 . . . and ^ |?/(6i)l = 1-855893 . . . have been used and M is given 
by eq. (|2.16p . The density in the degenerate limit vanishes where fi(r) does. 

Finally, for the velocity dispersion we find from eqs. (|2.14[) . (|2.16p and (|2.20[) 
The phase-space density at the origin Q(0) takes in the degenerate limit its maximal value 



4 

a deg (r) = <r de9 (0) ^(0 , a deg (0) = V O^fJ L || 



(2.23) 



In the quantum degenerate limit isq — > +oo, the DM fermions of the halo are in the minimally excited state, namely the 
(semiclassical) ground state for which the three basic quantities: halo radius, galaxy mass Mh and velocity dispersion, 
take their minimum values: 

/keV\^ / pc 3 Y / ^ km /keV\^ / pc 3 ^ 

rh,min = 24.51 . . . pc po — — , <T mi „(0) = 2.751 — p — — 

\ m J \ M Q J s \ m J \ M q 



.V/,,„,„ = 2.!):-!?)-! . . . I0 4 M Q ^-^j J Po |£ . (2.25) 



These minimum values are similar to the estimates for degenerate fermions eqs. (|1.13[) . as it must be. 

The masses of compact dwarf galaxies dominated by DM must be larger than the minimum mass Mh, m in eq. (12.25|) . 
The lightest known galaxy of this kind is Willman I (see Table H]). Imposing Mh,min < Mwuiman i = 2.9 10 4 M 
gives a minimal mass (a lower bound) for the WDM particle: 

m > m mm = 1.91 keV . (2.26) 

It must be stressed that these minimum physical values for degenerate WDM fermions eqs. (|2.25[) - (|2.26[) are indepen- 
dent of the WDM particle physics model. They are universal whatever is the shape of the distribution function f(E) 
in the non-degenerate regime. 

As can be seen from eq. (12.251) the minimal value of the WDM particle mass m min in eq. (|2.26[) is given by 

1 / V i 

/ 10 M G) \ I P n PC 

m min = 1.3094... keV ' ' ' 



I ^lightest galaxy J I fy[ 

As a consequence, the error on m m i n due to the observational error on the mass and central density of the lightest 
galaxy (Willman I, so far) is considerably reduced by the small exponents 1/4 and 1/8, respectively. 

X-ray galaxy observations provide upper bounds for the particle mass for a specific WDM candidate: the sterile 
neutrino. If the sterile neutrinos are described by the Dodelson-Widrow (DW) model this upper bound turns to be 

M 

m < 2.2 keV . (2.27) 

Therefore, for WDM formed by DW sterile neutrinos their particle mass must be in the narrow interval 1.91 keV < 
m < 2.2 keV. 

More precise data for ultracompact dwarf galaxies as Willman I will make our bound eq. l|2.26p more precise. 

Improvements on our lower bound eq. (|2.26[) as well as on the upper bound eq. (|2.27p from more precise galaxy data 
will lead to more stringent constraints on sterile neutrino WDM particle models. 
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B. Classical dilute limit: large galaxies 

In the classical dilute limit, v -C — 1, the FD distribution function eq. (|2.19[) becomes the Maxwell-Boltzmann 
distribution and we find from eqs. (12.121) . 

*d(2/ 2 -^) = e^ , ^(i/) = |V5Fe" , / 4 » = ^^e", (2.28) 



— 

(| V^J e" , i/(0) = , i/(0) - 



d 2 i/ 2 dv _ /3 A 

^(r) = ^=e^ =£. (2.29) 

In this limit from eq. (|2.28|) . the ratio lf(y)llf(y) becomes constant and the average velocity a(r) from eas. (|2.14p 
becomes uniform. Furthermore, the numerical resolution of eg. (12. 121) or eas. (|2.28l) in the classical t/n <C — 1 regime 
yields 

£ R = 37.11... e^ 3 > 1 and - i/(£r) = 0.05318 ... e Vo/3 < 1 , 

6> = 3.299... e"""/ 3 > 1 and - = 0.5578 ... e^ 3 < 1 . (2.30) 

Therefore, we see from eqs. (|2. 141) . (|2.16l ) that the size and mass for Uq -C — 1 and fixed po grow as 



M„ = 5.0202... e -™M o (!0^f| . r» = 61.71... e-/> (0 („ pc . 

We see that in the classical regime, vq — > — oo, and Mh can be arbitrarily large. Galaxies of any large mass can be 
obtained as solutions of the Thomas- Fermi equations (|2.12[) . 

C. Physical galaxy properties from compact dwarf galaxies to large galaxies 

The largest value for the phase space density corresponds to the quantum degenerate fermions limit while the 
smallest values appear in the classical dilute limit. 

In the left panel of fig. [5] we plot from eq. (|2.14p the dimensionless quantity 



& n(n , _ 7125 I|M /_m_* 
(keV)^ 3w a W 



In the right panel of fig. [2j we plot instead the dimensionless product 

£ =0.82296... 10»^iL . (2.32, 

M© V Po pc 3 4 (^o) V m / 

where M/j is the halo mass, namely the galaxy mass inside the core radius defined by eq. (!2.18[) and we used 
eqs. (|2.14j) - (l2.16[) . In both cases we consider the two values m = 1 and 2 keV and we put in the abscissa the product 

pc 3 V p c fkcVA 1 

— On = tin Eh 



rh y-j-j- poj = Ro 6i [~^~J in P arsecs ' ( 2 - 33 ) 

where r/, = Iq E,h is the core radius. The phase-space density Q(0) and the galaxy mass Mh are obtained 
by solving the Thomas-Fermi eqs. (|2.12p . We have also superimposed the observed values fi. 3 Qh/(keV) 4 and 
Mh\/MQ/[p a pc 3 ] (m/keV) 4 from Table [TJ Notice that the observed values Qh from the stars' velocity disper- 
sion are in fact upper bounds for the DM Qh- This may explain why the theoretical Thomas- Fermi curves in the left 
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FIG. 1: Density and velocity profiles, p(r)/po and cr(r)/a(0), as functions of r/lo for different values of the chemical potential 
at the origin vo- Large positive values of vq correspond to compact galaxies, negative values of uq correspond to the classical 
regime describing spiral and elliptical galaxies. All density profiles are cored. The sizes of the cores vh defined by eq. (|2.18[) are 
in agreement with the observations, from the compact galaxies where ru ~ 35 pc till the spiral and elliptical galaxies where 
th ~ .2 — 60 kpc. The larger and positive is vo, the smaller is the core. The minimal one arises in the degenerate case v» — > +oo 
(compact dwarf galaxies). 



pane of fig. [5] appear slightly below the observational data. Notice also that the error bars of the observational data 
are not reported here but they are at least about 10 — 20%. 

The phase space density decreases from its maximum value for the compact dwarf galaxies corresponding to the 
limit of degenerate fermions till its smallest value for large galaxies, spirals and ellipticals, corresponding to the 
classical dilute regime. On the contrary, the halo radius and the halo mass Mh monotonically increase from the 
quantum to the classical regime. 

Thus, the whole range of values of the chemical potential at the origin v§ from the extreme quantum (degenerate) 
limit > 1 to the classical (Boltzmann) dilute regime v$ <C — 1 yield all masses, sizes, phase space densities 
and velocities of galaxies from the ultra compact dwarfs till the larger spirals and elliptical in agreement with the 
observations (see Table [I]). 

From figs. [5] we can extract important information on the fermion particle WDM mass. 

We see from the left panel fig. [5]that decreasing the DM particle mass m moves the theoretical curves ft 3 Q(0) / (keV) 4 
towards smaller Q(0) values and larger galaxy sizes, one over each other. In the right panel of fig. [5] we see that 
decreasing the DM particle mass m displaces the theoretical curves (M^/M Q )^/ M@J[po pc 3 ] towards larger galaxy 
masses and sizes, one over each other. 

The small galaxy endpoint of the curves in figs. [2] corresponds to the degenerate fermion limit vq — > +oo and its 
value depends on the WDM particle value to. For increasing to, the small galaxy endpoint moves towards smaller 
sizes while for decreasing to, it moves towards larger sizes. 

We see from figs. [5] that decreasing the particle mass beyond a given value, namely for particle masses to < 1 
keV, the theoretical curves do not reach the more compact galaxy data. Therefore, to < 1 keV is ruled out as WDM 
particle mass, in agreement with the bound eq. (|2.26[) . 

For growing to > keV the left part of the theoretical curves corresponding to the lower galaxy masses and sizes, will 
not have observed galaxy counterpart. Namely, increasing m ^> keV would show an overabundance of small galaxies 
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(p pc 3 /M ) 1/6 r h [pc] (po pc 3 /M Q ) 1/6 r h [pc] 

FIG. 2: In the left panel we display the galaxy phase-space density ft 3 Q(0)/(keV) 4 = (m/keV) 4 (VT25/3 tt 2 ) [if (v )] 
defined in eq. (|2.14| obtained from the numerical resolution of the Thomas- Fermi eos. (|2.12p for WDM fermions of mass m = 1 

o 1 4 

and 2 keV versus the ordinary logarithm of the product log 10 {rh [pc po/Mq]«} = log 10 [i?o £,h (kev/m)^] in parsecs. The 
red stars * are the observed values of h 3 Q(0)/(keV) 4 from Table [T] Notice that the observed values Qh from the stars' 
velocity dispersion are in fact upper bounds for the DM Qh and therefore the theoretical curve is slightly below them. In the 
right panel we display the galaxy mass (M/M@ ) \J Mq / [p p pc 3 ]" = 0.82296 10 5 (keV/m) 4 & \v'(£h)\ /-f^Vo) obtained from 
the numerical resolution of the Thomas- Fermi eqs. (|2.12[) for WDM fermions of mass m = 1 and 2 keV versus the product 

rh [pc 3 po/Mq}~s — R £ h (keV/m) 3 in parsecs. The red stars * are the observed values of (M/Mq)^/Mq/[po pc 3 ] from Table 
HI Notice that the error bars of the observational data are not reported here but they are at least about 10 — 20%. 



(small scale structures) which do not have counterpart in the data. This is a further indication that the WDM particle 
mass is approximately around 2 keV in agreement with earlier estimations [ToL [l2| . 

In addition, the galaxy velocity dispersions from eq. (|2.14l) - (|2.16l) turn to be fully consistent with the galaxy obser- 
vations in Table |U 

The semiclassical Thomas-Fermi approach provides much stronger results than just the quantum bound on the 
phase-space density Q(r) only based on the Pauli principle, 

4 

Til 

Q{?)<K^ , K = 0{1). 

The primordial phase-space density Q fulfils this quantum bound Through classical time evolution, Q(r) in 
average (Q) will always fulfil this quantum bound since as it is known, Q can only decrease with time. This quantum 
bound gives lower bounds in the galaxy core sizes in the range ~ 0.1 pc 9] which are of the same order of magnitude 
as the core sizes obtained in the classical (non-quantum) A^-body WDM simulations [l5j], which are unrealistically 
small cores compared with the observations. On the contrary, the Thomas-Fermi approach includes the quantum 
pressure and indeed succeeds to provide galaxy cores with the right sizes as shown in this paper. 

To conclude, cqs. (|2.14| - (|2.16p indicate that the galaxy magnitudes: halo radius, galaxy masses and velocity dis- 
persion obtained from the Thomas-Fermi quantum treatment for WDM fcrmion masses in the keV scale are fully 
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consistent with all the observations for all types of galaxies (see Table |T| . Namely, fermionic WDM treated quantum 
mechanically (as it must be) is able to reproduce the observed sizes of the DM cores of galaxies. These results strenght 
the discussion in sec. II B I that compact galaxies are supported against gravity by the fermionic WDM quantum pres- 
sure. 

It is highly remarkably that in the context of fermionic WDM, the simple stationary quantum description provided 
by the Thomas-Fermi approach is able to reproduce such broad variety of galaxies. 

In addition, WDM simulations produce the right DM structures in agreement with observations for scales > kpc 

a- 
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